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1 Calculating Absorption of a Lossy Gas Using Q 


Calculating microwave opacity from a weakly absorbing gas mixture using a resonator requires 
measuring the quality factor of that resonator which necessitates accurately determining the center 
frequency (/ 0 ) and the half power bandwidth (A/) of a noisy resonant line. The center frequency 
can be determined very accurately and varies very little over many measurements (a few kHz at 
GH* frequencies or a few hundredths of a percent). The greater source of error in estimating the Q 
of a resonator comes from the bandwidth measurements. The half power bandwidth is determined 
essentially by eye-fitting a curve over a noisy resonant line and measuring with a spectrum analyser. 

Assuming uncorrelated measurement error (an unbiased observer) the best estimation of the 
bandwidth is obviously the mean of many measurements. Therefore, 

(1) 

4=1 

where A/, are the individual bandwidth measurements. To determine the accuracy of the measure- 
ments we calculate the sample variance, 
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and finally the variance of our estimation from the “true* bandwidth due to electrical noise, 
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where is the “student-t* for the 1 a confidence level given a Gaussian distribution. Recall that 
the student-t is a distribution to characterise the confidence level of a finite sample set where the 
degrees of freedom are the number of samples. For ten measurements, t a = 1.1 (5, pp. 255, 260]. 

One can also calculate some figures of merit to determine the stationarity of the data [5, pp. 
255-6]. Divide the ten measurements into two groups of five and calculate the sample variance 
(s^Sj) and mean (mi, m 2 ) of each group where and mi > mj. We can define a figure of 

merit for stationarity of the variance as F 0 = 0.156sJ/s^. If Fp< 1 t hen the variance is likely to 
be stationary. For the mean, we define t Q = 0.974(mi - ma)/>/a \ + s\ M figure of merit and if 
t G < 1 then the mean is likely to be stationary. If either one of these conditions does not hold the 
measurement is not stationary and should be considered suspect. 

There are also instrument errors which must be taken into consideration as welL These stem 
from the limited accuracy of the spectrum analyser and are calculated as follows [2]: 

o 0 w 10 _6 /o + 0.15RBW + 0.05SPAN + L c + 250 Hi (4) 


o 

< 

< 

O 

c 

tO 

r v 

cn 

L. 

o o 

CO Q; 
rn O 
m w 

0 

t z 

CL O 
O i-» 

1 K 
< u 
to 

<o a 

-Z U j 

^ OC o 



B2 


a A « 10‘ 6 BW + 200 iV + 4 Lew Hz (5) 

where a Q i* the standard deviation of the center frequency measurement due to the spectrum ana- 
lyser accuracy, a A the standard deviation of the bandwidth measurement, RBW is the resolution 
bandwidth, L the least significant digit, BW the bandwidth and N the harmonic number (see tables 
1 and 2). 
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RBW |Hs| 

300 

1000 

30000 

SPAN [kHz 

50 

100 

2000 

N 

1 

2 

4 

Lo 

10 

100 

1000 

Lbw 

10 

100 

1000 


Table 1: Spectrum analyser parameters 



o 0 [kHz] 

(7 A [Hz] 
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25 

240 

X 

90 

800 

K 

320 

4815 


Table 2: Instrument uncertainty for the frequency bands. 


It is then possible to estimate the measured bandwidth and characterize its uncertainty. What we 
measure from a conventional spectrum analyzer, however, is not the true spectrum of the resonator 
but rather the convolution of the resonator spectrum and the sweep filter spectrum. If Y(f) is the 
output, X(f) the input and H(f) the filter, then 

Y(f) = H(f)*X(f). (6) 

In the time domain, we have y(t) = h(t)x[t) where y(t) is the inverse Fourier transform of Y(f) and 
so on. Therefore, the signal spectrum is, 

*(/) = W)/M01 ( 7 ) 


where T[-\ denotes a Fourier Transform. If we assume a Gaussian filter and a Gaussian input then 
the actual signal bandwidth is simply 


A/q = \/a/ 2 -RBW 2 


( 8 ) 


where A/q is the “true* half power 
of the resonator is then simply 


bandwidth and Af is the measured mean bandwidth. The Q 
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We will discuss the error in the next section. 
To calculate the absorption recall that 
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( 10 ) 
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Figure 1: a/ fi vs. e"/e' for full expression (eq. 10) and approximation (eq. 11). 


where $ = 2r/X (with the gas present), a is the absorption and e = e'-je" the complex permittivity. 
Since c" C e* for gases in the microwave region this reduces to 


a e" 

— fa — . 
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Note that this is field attenuation — to get power attenuation one must multiply this by 2. 
The Q of a dielectric, Q d) is given by 
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where the approximation holds for a gas at microwave frequencies. Therefore, we see that 

fi * 1 
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2 Q d *Qd 


To determine Q d we must examine the Q of a lossy, coupled resonator. 

The measured quality factor of an evacuated resonator can be given by |l, p. 404] 
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and in the presence of a lossy gas 
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where 
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(15) 
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QTg Qcg Qd Qeg(l) Qeg( 2 ) 

Qrv.g — the total (measured) Q without and with the gas 
Qcv.g — cavity Q without and with the gas 
Qev.g(i,2) — external Q on both sides without and with the gas 
Qd = Q of the gas. 

Since a symmetrical cavity is used, assume that Q es ,„(i) = Qtg,v{ 2) = Qeg.v Additionally [ 1 , p. 
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where t g and t v are the transmissivities of the resonator at resonance with and without the gas 
respectively (assuming a lossless gas). Substituting into equations 14 and 15, subtracting and solving 
for 1/Qd yields 

_L _ l ~ >/*g _ * ~ L_ (17) 

Qd Qtq Qtv Qcv Qcg 

Note that Qcv,g takes conductivity losses and diffraction losses into account. The Q for finite 
conductivity, Q a> of a circular cavity varies as l/y/J [3, p, 352] we see therefore that 


«,/4wi. 


(18) 


Q 


For diffraction, likewise, the ratios of the Q’s will be on the order of the ratios of the wavelengths 
and also equal to about 1, therefore Qcv ** Qcg • Substituting these approximations into 17 yields 


(19) 


_L _ 1 ~ y/^ _ 1 ~ yfc 
Qd Qt 9 Qtv 

Recall that t g assumes a lossless gas. Actually, the measured transmissivity, t g fJlra ,, equab t g c~ ad ' lf 


or 
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( 20 ) 

where d c jj = QA/2* is the effective path length of the resonator. We can then solve iteratively for 
the actual a. 

We therefore see that 




where A is in meters. In dB, we can express a as 


aaa = 8.686^ «/»■ 


( 21 ) 


( 22 ) 


The above analysis neglects the effects of dielectric loading [4]. Using an amount of lossless 
gas with the same refractivity value as the lossy gas in lieu of merely an evacuated resonator will 
remove the effects of the dielectric loading, provided you are far enough from the line center of the 
absorption line so that anomalous dispersion is not present. 


2 Calculating the Error 

Define a quantity 


t = 9(fog, A /„ /«,, A/.) = ^ 
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where A f g and A /„ are the loaded and unloaded half power bandwidths and and f ov are the 
loaded and unloaded center frequencies. We define this quantity since we wish to ignore the 1/A 
dependence in front of a in the interest of symmetry. Then to first order 
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Note that this ignores variation in 7 = 1 — y/t. 

Assume that < 6£ >— 0 that is, we are dealing with sero-mean processes; therefore 

a \ =< s? >=< r; > + < v] > -2 < r,r„ > (25) 

where 
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The bandwidth variation consists of instrument accuracy and noise; i.e. 6 A fa = SA/sa + SAf^i, 
Throughout, we will assume that electrical noise and the spectrum analyser accuracy are uncorre- 
lated. Since we have assumed sero-mean processes and uncorrelated accuracy /noise then 


<Sfl> = (29) 

<6Aft> = o\+o 2 Ni . (30) 

The uncertainty of the measured gas absorption, a, is then 

± dB/m (31) 

where we have neglected the uncertainty in the measurement of A. There are three cases that we 
will consider: the uncorrelated case, the “worst" correlation case and the “best* correlation case. 


2**1 Uncorrelated Case 


In the uncorrelated case, where naturally all variations are uncorrelated with one another, we find 


that 
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2..2 Worst Case 

In the worst case scenario (the greatest error), 6 f v 6Af g and 6f g 6Af v terms are completely correlated 
while the other terms are completely anti-correlated and, of course, the electrical noise terms are 
uncorrelated with the instrument accuracy terms. This yields 
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This is the case typically used for the error, since a few small error terms have been neglected in 
this derivation. 



2. .3 Best Case 

The beet case scenario is the reverse of the worst case; that is 6f v 6&f g and 6f g 6&f v terms are 
completely anti-correlated while the other terms are completely correlated. This yields 


< r? > 
< r„r„ > 


+ <r| + ffj/f - 


2(7</rA ] 

Qi \ 

OqO& VqG A 1 

Q„ Qv \ ' 


(36) 

(37) 


References 

[1] Ginzton, E.L., Microwave Measurements, 1957, McGraw-Hill Book Co., Inc. 

[2] HP Manual for HP8562B 

|3] Posar, D. M., Microwave Engineering, 1990, Addison- Wesley Publishing Company, Inc. 

|4] Spilker, T.R., New Laboratory Measurements on Ammonia’s Inversion Spectrum, With Impli- 
cations for Planetary Atmosphere’s, J. Geophys. Res., 1993, 98:5539 

[5] Warner, F.L., Microwave Attenuation Measurements, 1977, Peter Peregrinus Ltd. 



